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Abstract 

The  classical  theory  of  longitudinally  slot- 
ted vails,  which  substitutes  an  approximate  homo- 
geneous vail  boundary  condition  for  the  true 
siixcd  conditions,  is  extended  in  several  rcsjxicts. 
Based  on  recent  experimental  findings  at  the  FFA, 
an  inviscid  flow  model  is  adopted  in  which  the 
outgoing  slot  flow  penetrates  into  the  plenum 
chamber  as  a thin  Jet,  wfiile  the  re-entering 
flow,  admitting  quiescent  air  from  the  plenum 
chamber  into  the  test  section,  induces  a longi- 
tudinal separation  bubble  at  plenum  pressure 
along  the  slot  and  adjacent  parts  of  the  test 
section  wall.  The  three-dimensional  analysis, 
based  on  the  assumption  that  the  slots  are  narrow, 
retains  quadratic  cross-flow  terms  in  tho  pres- 
sure equation  and  allows  the  slots  to  be  few  in 
number  and  have  non-uniform  distribution  and  ge- 
ometry. A family  of  homogeneous  boundary  condi- 
tions is  obtained,  each  of  successively  higher 
a' curacy.  Application  to  tho  design  of  inter- 
forcncc-freo  transonic  test  sections  is  discuss- 
ed. Unsteady  effects  are  also  considered. 

1.  Introduction 

It  is  a remarkable  fact  tliat  theory  plays  a 
very  minor  roll  in  the  design  and  use  of  slotted 
transonic  test  sections,  in  particular  when  one 
considers  that  the  theory  was  already  well  devel- 
oped at  an  early  stage1'3.  There  might  be  sewjnil 
reasons  for  this  but  one  of  them  seems  to  be  that 
the  flow  models,  mainly  inviscid  ones,  had  not 
boon  tested  by  careful  experiments,  while  incon- 
sistencies were  alleged  to  appeal'  in  trying  to 
apply  the  theoretical  results  to  practical  wind 
tunnel  flows.  Therefore,  it  seems,  reliance  was 
placed  on  empirical  methods,  determining  the 
slot  width  so  as  to  minimize  choking  effects 
around  Mach  one  and  accepting  as  free-stream 
conditions  the  upstream  conditions  in  the  empty 
test  section  as  calibrated  against  the  plenum 
pressure.  It  became  a widespread  belief  that 
viscous  effects  of  an  unknown  and  complicated 
nature  are  present  in  the  slot  flow  and  tliat 
therefore  on  empirical  approach  is  all  there  is 
available.  In  consequence  the  development  of 
the  theory  came  to  a virtual  standstill. 

In  a recent  experimental  investigation  of 
•lot  flows4  it  has  been  found  tliat  under  typical 
test  conditions  tho  slot  flow  is  not  necessarily 
dominated  by  viscosity  and  that  therefore  one 
can  define  an  inviscid  flow  model  which  is  rea- 
sonably closo  to  reality.  This  model,  which  is 
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not  as  simple  as  aoaie  of  the  early  ones,  will  be 
set  down  in  the  next  section.  For  a fuller  dis- 
cussion of  the  problems  involved  see  Ref.  4 and 
references  given  therein. 

Our  present  task  is  to  base  on  this  model  a 
general  inviscid  theory  for  the  wall  interfer- 
ence caused  by  a slotted  wall.  In  the  classi- 
cal theory  of  this  kind5*3,  based  on  the  slonder- 
body  approximation,  a description  is  obtained  of 
the  combined  effect  of  the  mixed  boundary  condi- 
tions at  the  slots  and  slats  in  the  form  of  a 
far  simpler,  homogeneous  boundary  condition,  re- 
lating the  average  pressure  difference  across 
the  wall  to  the  average  streamline  curvature 
normal  to  the  wall.  The  simplification  is  made 
possible  by  postulating  a large  number  of  simi- 
lar slots  uniformly  distributed  over  the  inter- 
fering wall. 

In  order  to  minimize  any  viscous  effects, 
and  possibly  for  other  reasons,  it  is  desirable 
to  keep  the  number  of  slots  small.  Furthermore, 
if  the  slots  are  few  in  number  their  proper  lo- 
cation might  become  important  when  minimizing 
the  wall  interference.  Consequently  there  is  a 
great  need  for  freeing  the  classical  theory  from 
its  inherent  restrictions  while  keeping  its  sim- 
plicity. Wo  sliall  achieve  this  by  employing  a 
modified  method  of  approximation,  based  on  the 
much  less  restrictive  assumption  that  the  slot 
width  is  sna’  F^ompared  to  the  distance  between 
slots.  From  the  case  of  uniform  slot  distribu- 
tion tills  is  known  to  be  a workable  approach*'10. 
In  the  present  more  general  context  it  leads  to 
a straightforward  application  of  the  method  of 
■latched  asymptotic  expansions. 

In  our  Analysis  we  shall  aim  at  a fair  level 
of  generality.  This  might  for  example  facili- 
tate later  inclusion  of  corrections  for  viscous 
effects.  The  main  concern  will  be  with  three- 
dimensional  tests,  the  case  of  two-dimensional 
tests  being  already  treated4.  Also,  we  shall 
not  excludo  the  case  of  unsteady  flow,  which  is 
clearly  quite  important  in  connection  with  os- 
cillatory testing11 . However,  in  the  present 
paper  the  applications  of  the  theory  will  be  re- 
stricted to  a few  very  simple  cases. 

Tho  theory  turns  out  to  produce  not  one  but 
a whole  family  of  possible  homogeneous  boundary 
conditions,  each  corresponding  to  a specific 
degree  of  resolution  of  the  details  of  tli*  wall 
flow  as  ’filtered*  through  the  averaging  proce- 
dure. Tho  choice  for  any  particular  application 
will  have  to  be  based  on  considerations  of  ac- 
curacy. There  is  also  some  freedom  In  the  me- 
thod of  applying  the  boundary  conditions,  de- 
pending on  the  kind  of  problem  to  be  solved.  If 
the  slots  are  to  be  adjusted  for  zero  interfer- 
. once  by  employing  measured  wall  pressure  distri- 
butions, a somewhat  different  type  of  boundary 
condition  will  ho  required  from  if  the  adjustment 
is  to  be  based  on  a p re -computed  interfor-'iice- 
J free  flow  field.  We  shall  return  to  this  point 
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ng.  1.  Possible  inviscid  flow  patterns  (not  to  scale). 
OH 8®  = fast  air  from  the  slot. 


in  later  soctions  as  well  as  to  the  i elated 
questions  of  how  to  define  and  compute  wall  inter* 
Terence  corrections  and  how  to  use  test  section 
calibrations. 


In  order  to  translate  the  slot  flow  model  in- 
to a set  of  boundary  conditions  to  be  used  with 
the  flow  equations  further  simplifications  must 
be  introduced.  First  of  all  there  is  no  need  to 
consider  in  detail  the  development  of  the  jet  in- 
side the  plenum  chamber.  From  the  point  of  view 
of  the  slot  flow  it  should  be  sufficient  to  des- 
ignate a surface  Spo  across  the  slot  exit  on 
which  the  plenum  pressure  may  be  taken  to  act. 

It  must  also  be  specified  how  the  fast  air  in 
the  slot  is  split  off  to  return  to  the  test  sec- 
tion: making  the  obvious  choice,  we  shall  as- 

sume tliat  the  downstream  end  of  Sp0  coincides 
with  the  upstream  end  of  the  free  surface  S-  , 
which  is  the  boundary  between  the  fast  air  and 
the  quiescent  plenum  air  (see  Fig.  l).  Tlius 
and  Sp  together  form  a surface  on  which  we  have 
plenum  pressure.  Obviously  thore  is  no  need,  nor 
any  real  possibility,  to  dotennine  Spo,  and  hqioo 
Sp  in  its  subsequent  development,  with  any  ac- 
curacy. We  shall  interpret  this  as  a licence  to 
stake  a choice  which  renders  the  analysis  simple. 


2.  An  Inviscid  Model  for  the  Slot  Flow 


The  slot  flow  etodel  to  be  adopted  is  present- 
ed in  Fig.  1.  In  the  upstream  part  of  the  test 
section  the  flow  is  going  outwards  through  the 
slot  into  the  plenum  chamber,  where  it  forms  a 
thin  jet.  The  flow  inside  the  slot  is  attached, 
as  Indicated  in  cross  section  (a),  and  separa- 
tion occurs  at  the  sharp  slot  edges  on  the  pie  nun 
side.  Above  the  model  the  slot  flow  turns  back, 
leaving  the  jet  to  continue  on  its  own  into  the 
interior  of  the  plenum  cliambcr.  This  'splitting' 
of  the  fast  air  into  two  separate  streams  is 
shown  beginning  at  section  (b).  At  (c)  the  fast 
air  in  the  slot,  having  entered  farther  upstream, 
is  returning  to  the  test  section.  At  (d)  the 
fast  air  has  left  the  slot  and  behind  it  appears 
a 'bubble'  of  quiescent  air  at  plenum  pressure, 
the  boundary  of  which  is  expanding  into  the  test 
section  flow.  Typically,  the  bubble  is  narrow 
and  extends  along  the  slot.  Presumably  it  col- 
lapses onto  the  slot  farther  downstream  if  the 
cross  flow  turns  back  again  towards  the  wall. 

What  happens  if  it  is  struck  by  a shock  wave  from 
the  model  is  not  known. 


Having  so  far  tacitly  assumed  that  the  flow 
is  steady,  we  must  also  consider  how  the  flow 
model  can  be  generalized  to  become  applicable  to 
unsteady  flows.  Obviously,  the  free  surface  S 
Bust  be  allowed  to  move.  The  possibility  that 
pressure  waves  propagate  inside  the  plenum  clia m— 
ber  Bust  also  bo  considered.  Consequently,  the 
pressure  cn  SpQ  + Sp  cannot  be  takon  to  be  known 
in  advance.  Trie  so  are  serious  complications 
which  can  perhaps  not  be  handled  without  making 
further  simplifications,  sucli  as  assuming  the 
unsteadiness  to  be  a small  porturbaticn  of  a 
steady  flow. 


This  is  tho  picture  of  the  slot  flow  arrived 
at  in  Ref,  4.  The  description  leaves  undecided 
whether  the  high  speed  air  returning  from  tho 
•lot  to  tlio  test  section  is  a vorticity-carrying 
■lug,  as  in  (cl)  and  (dl),  or  whether  it  expands 
around  the  slot  edges  without  separation,  as  in 
(c2)  and  (d2).  The  experimental  evidence  in  Ref. 

although  not  quite  conclusive,  points  to  the 
fonsor  alternative.  It  might  well  be  that  both 
types  of  flow  may  occur.  In  order  to  avoid  com- 
plications at  this  slago  wo  shall  assume  the 
second  type  of  flow,  but  when  developing  the 
analysis  wo  shall  keep  the  alternative  in  mind 
as  well  as  tlio  possiblo  need  for  viscous  correc- 
tions. 


itions 


The  test  section  wall,  before  the  slots  are 
introduced,  is  taken  to  be  a cylinder,  Sw 
(Fig.  2).  Tlio  longitudinal  slots,  numbored 
1 to  N,  are  connected  to  plenum  chambers  with 
quiescent  air  at  prescribed  pressures  pJD(i  ■ 
1,2,..,N)<  Tho  x-uxia,  parallel  to  S , points 


In  tho  flow  direction.  Tho  radius  vector  r Is 
orthogonal  to  the  x-axis.  Tim  hydraulic  rrulius 
of  Sw  Is  employed  as  the  unit  of  length. 


Wg.  2.  Slotted  test  section. 


The  model  and  its  wake  ore  located  within 
'S  the  cylinder  S„,  , parallel  to  Sw  . The  flow  is 
1_  taken  to  be  inviscid  between  Sm  aiid_  , 'as 
'I)  well  as  in  the  slots  and  plenum  chomtjors.  The 
jl  flow  of  fast  air  will  be  described  by  a small 
• irrotatioual  perturbation  of  a uniform  reference 
flow  parallel  to  the  x-axis  (at  density 
pressure  , velocity  U and  Mach  number  M), 
tlie  perturbation  to  bo  considered  as  produced  by 
tho  distribution  of  normal  velocity  set  up  by  the 
model  on  Sm  . This  distribution  is  of  course 
not  known  in  advance.  Our  ignorance  of  the  flow 
Inside  Sm  is  the  reason  for  the  wind  tunnel 
test  to  begin  with,  but  from  the  shape  of  the 
model  one  can  often  estimate  it.  with  suffi- 
cient accuracy  for  computing  the  wall  interfer- 
ence. This  we  assume  to  be  so  and  choose  for 
the  reference  state  the  free-stream  state  at 
which  the  estimate  is  mode. 

, Imt  U*  <p(x,r,t)  be  the  perturbation  veloc- 
ity potential,  normalized  to  be  zero  in  tho  ref- 
erence flow.  To  first  order  for  transonic  flow 
It  satisfies  the  following  differential  equation 
between  Sj,  and  the  outer  boundary  of  tlm  fast 
air  i 

= [M,-1+Ma(Y+l)<Px]«Pxx  + ^'(>xt+^Ttt.  0) 

Here  A is  the  Laplacian  in  planes  x = constant, 
while  T is  the  ratio  of  specific  heats  in  the 
reference  state.  The  inner  boundary  condition 
for  9 is  at  a point  P 

«*  s.  « 1n  *=  , (2) 

where  n denotes  differentiation  In  the  normal 
direction  and  F is  the  normal  velocity  distri- 
bution (which  we  have  assumed  known).  The  outer 
boundary  condition  on  solid  surfaces  adjacent  to 
fast  air  ia  similar: 

on  S between  slots:  ■ «•  0 , 

w ’n 

Internally  in  slots:  9n«  H(p)  . (3) 

H vanishes  aherever  the  walls  of  the  slots  are 
parallel.  (More  generally,  we  could  easily  al- 
*Pn  to  be  non-vanlshlng  on  Sw  as  well, 
thus  accounting  for  small  deviations  from  cylin- 
drical geometry  as  well  as  for  tho  displacement 
effect  of  wall  boundary  layers.) 

The  remaining  outer  boundary  condition  in- 
volves the  pressure.  In  the  present  approxima- 
tion the  Bernoulli  equation  ia 


P-R»  * “ Q*,  l^(jj  9t«-9x*  5 Grad*  9 ) , («•) 

where  Grad  denotes  the  gradient  operator  in  a 
plane  x « constant.  Tho  quadratic  cross-flow 
term  is  needed  only  In  the  slot  regions,  whore 
tho  cross-flow  velocity  might  be  considerably 
larger  than  elsewhere.  In  the  case  of  steady 
flow,  the  jet  and  free-surface  boundary  con- 
dition therefore  takes  the  form 

on  Sp<>  ♦ Sp  : Grad*?  = - 6 U)  , 

p ul-  p 

•>t)c  P~— ~ (i  = 1.....N)  ♦ (5) 

The  corresponding  condition  for  unsteady  flow 
will  bo  discussed  later. 

In  order  to  formulate  an  upstream  boundary 
condition  we  shall  assume  that  tho  slots  alL  be- 
gin at  x = xQ  and  that  their  widths  increase 
smoothly  from  zero.  With  a stretch  of  parallel 
walls  between  the  contraction  and  the  beginning 
of  the  slots  — since  F is  likely  to  be  small 
that  far  upstream  — the  flow  at  x = xQ  ought 
to  be  uniform  (except  for  receding  waves  in  the 
unsteady  case).  We  shall  in  fact  take  this  to 
bo  so  and  prescribe  as  an  upstream  boundary  con- 
dition for  steady  flow  that  <p  takes  a constant 
valuo,  tp0  say,  on  the  plane  x = At  the  same 

time  wo  must  adjust  F to  vanish  at  xQ.  Note 
that  the  Mach  number  M0  of  the  entrance  flow 
is  a controllable  test  parameter  which  deter- 
mines a particular  valuo  for  <px  . We  shall  not 
specify  any  downstream  boundary  condition,  only 
assume  that  the  slotted  part  of  the  test  section 
is  long  enough  for  the  conditions  at  the  down- 
stream end  not  to  influence  the  flow  at  the 
model. 

In  the  present  approach  the  wall  interfer- 
ence is  obtained,  obviously,  by  subtracting  from 
9 on  Sm  the  corresponding  distribution  ob- 
tained with  unbounded  flow  outside  Sm  (using 
the  same  reference  flow)  in  precisely  the  situa- 
tion prescribed  when  estimating  F . If  we  are 
to  minimize  the  interference  we  must  explore  the 
influence  on  9 of  the  several  test  parameters 
at  our  disposal.  The  mixed  form  of  the  set  of 
outer  boundary  conditions  for  9 constitutes  a 
major  difficulty  when  performing  this  task. 

fr.  The  Method  of  Approximation 

In  order  to  overcome  this  obstacle  we  shall 
introduce  an  approximation  <p(x,r, t)  for  <p  , 
satisfying  the  same  dif fcrential~ equation  and 
the  same  inner  boundary  condition,  but  a new 
outer  boundary  condition.  This  boundary  condi- 
tion ought  to  be  as  simple  as  possible  from  tho 
point  of  view  of  computing  9 , consistent  with  tie 
requirement  that  9 must  be  closely  equal  to  9 
on  Sm  (whore  the  interference  is  to  be  computed). 
Thus,  the  new  boundary  condition  shall  be  re- 
quired to  be  homogeneous  and  local  in  the  sense 
that  it  is  a regular  functional  relationship 
between  9 and  its  normal  derivative  over  the 
entire  boundary  Sw,  while  values  of  9 or  9 
insido  Sw  , e.g.  on  Sm  , must  not  bo  explicit- 
ly present.  Clearly  the  boundary  conditions  for 
9 are  not  homogeneous  in  this  sense,  although 
they  arc  local. 


The  potentials  5 and  tp  nro  expected  to 
be  nearly  equal  almost  evei-yvhere.  In  particular 
In  the  neighbourhood  of  Sm  , and  to  be  essential- 
ly different  only  wticre  slots  aro  located.  Tliere 
the  cross-flow  derivatives,  hut  not  the  poten- 
tials themselves  or  their  derivatives  with  re- 
spect to  x or  t , aro  expected  to  be  much  dif- 
ferent, Noting  ttiat  tlie  differential  equation 
(l)  contains  cross-flow  derivatives  only  In  the 
left-hand  memlier,  wo  arrive  at  our  basic  method 
of  approximation!  we  neglect  the  difference 
between  $ and  9 in  the  night-hand  member  and 
postulate  that 

1 (»-*)■  0 . (6) 

This  Is  of  course  the  alcnder-body  approximation 
applied  to  the  slot  flow.  There  is  additional 
support  for  its  validity  in  that  in  transonic 
flow,  and  at  low  frequencies,  all  terms  in  the 
right-hand  member  of  ( 1 ) are  small.  Hut  this 
also  points  to  the  danger  tliat  the  approximation 
(6)  may  not  bo  very  good  near  shock  waves  and 
high-f roquoncy  receding  waves.  At  best  the  ap- 
proximation can  be  verified  a posteriori. 

The  differential  equation  (6)  can  be  inte- 
grated immediately  to  give 

• 9«9«-0-0,  (7) 

where  0(rjx,t)  and  5(rjx,t)  are  two-dimen- 
sional harmonic  functions  satisfying  the  santo 
type  of  boundary  condition  as  9 and  9 . More 
specifically  0 is  taken  to  satisfy  an  inner 
condition  (2)  with  the  normal  derivative  pre- 
scribed as  f(P,t)  on  S,„  (f  being  similar 
to, but  in  general  different  from  F)  and  in  addi- 
tion to  satisfy  the  condition  (3)  for  the  normal 
derivative  at  the  outer  boundary  as  well  as  the 
conditions  (5)  for  the  slot  pressures  (or  a cor- 
responding sot  for  unsteady  flow).  In  general 
this  will  determine  0 uniquely  in  terms  of  f • 
Similarly,  0 is  taken  to  satisfy  the  same  inner 
boundary  condition  for  the  normal  derivative  as 
0 and  a new  homogeneous  and  local  outer  bound- 
ary condition  chosen  so  as  to  make  0 uniquely 
determined  and  easy  to  calculate,  rendering  it 
at  the  same  time  closely  equal  to  0 on  Sm  . 

If  the  outer  boundary  condi  tion  for  0 1ms 
the  required  property  of  not  containing  f ex- 
plicitly, then  it_might  be  eligible  as  the  con- 
dition defining  tp  . Let  us  assume  this  to  be  so 
and  solve  Eq.  (l)  for  9,  applying  the  inner 
boundary  condition  (2)  together  with  the  outer 
boundary  condition  thus  taken  over  from  0 . 

This  caui  bo  done  without  specifying  f . We  then 
conclude  from  (7),  assuming  the  basic  approxima- 
tion (6)  to  be  true,  that  9 is  closely  equal 
to  <p  on  Sm  , as  it  should  bo. 

The  crucinl  question  is  now  whether  9,  ns 
approximated  by  (7),  satisfies  the  outer  boundary 
conditions  (3)  and  (5)  to  sufficient  accuracy. 

To  verify  this,  choose  f so  as  to  make  0|t  = ^ 
on  S„  . The  equality  or  0 and  9 on  Sw  (in 
virtue  of  the  common  boundary  condition)  is  there- 
by extended, • approximately,  to  a neighbourhood 
of  Sw  • Therefore,  in  consequence  of  (7),  it 
might  be  assumed  that  in  n similar  neighbourhood 
(which  is  taken  to  include  the  slots)  9 ~ 0 , 
ffx  - tK  , «Pt  — 0t  * Grad  9 — Grad  0 , and  also 
. that  the  free  boundaries  Sp  coincide.  Then 


9 satisfies  the  boundary  conditions  (3)  and  (5), 
and  the  problem  of  computing  9 on  Sg,  , as  in- 
fluenced by  tho  test  section  wall,  has  lieen  re- 
duced to  the  simpler  problem  of  computing  9 , 

Note  that  in  this  process  it  might  be  un- 
necessary to  compute  0 , f , or  $ , since  none 
of  them  is  present  in  the  boundary  condition  for 
9.  However,  If  we  want  to  determine  conditions 
at  tho  wall,  for  example  wall  pressures,  then  we 
shall  need  to  know  0 there. 

5»  Conformal  Mapping  and  Analytic 
Continuation.  Smoothing. 

In  order  to  arrive  at  tne  outer  boundary  con- 
dition for  9 we  must  analyse  0 and  0 . They 
are  both  two-dimensional  harmonic  functions  in 
regions  with  shapes  independent  of  x . It  is 
therefore  natural  to  perform  a conformal  mapping 
of  the  annular  region  between  Sw  and  Sm  in 
each  plane  x = constant,  mapping  the  boundaries 
onto  concentric  circles.  In  addition  to  other 
simplifications  this  will  permit  us  to  reform- 
ulate the  inner  boundary  conditions  for  0 and 
? by  analytic  continuation  to  bo  imposed  at  the 
centre  of  the  circles. 

The  mapping  is  scaled  to  leave  the  cross- 


sectional  area  of 


invarian  t . 


be  polar  coordinates  in  the  transformed  plane. 
Hence  Sw  is  mapped  onto  the  unit  circle. 

We  now  postulate  that  in  analysing  $(r,9;x,t) 
it  is  sufficient  to  include  only  terms  up  to  and 
including  the  order  u in  a Fourier  expansion 
'with  respect  to  0 . This  is  a decisive  step:  it 
specifies  a 'filler'  which  permits  us,  by  choice 
of  v , to  approximate  0 by  $ with  control] a— 
bio  smoothness  at  the  wall  and  precision  at  the 
model.  Then  0 must  have  the  form 

$ = A In  r+B  +£  T (A  cosjQ  ♦ B .sin jO  )r~^  + 

J=1L  J J 

♦ (DjCosj9+ EjSinj9)rJ]  , (8) 

where  the  coefficients  A , A,,  B.,  S and  E 
are  all  functions  of  x and  ^t.  J 

The  corresponding  expression  for  0(r,9;x,t), 
having  the  same  singular  structure  at  tho  origin, 
is 

0 = A in  r+D  +£  [ (A  .cos  j9  ♦ B.sinj9)r-^  «• 

. 0 °j=tL  J J 

♦ (DjCosj9  + EjSinJ9)r^J  ♦ 6 , (9) 

where  tho  remainder  S , containing  the  harmonic 
components  required  for  describing  details  of 
the  slot  flow,  is  0(ru+  ) as  r-*0.  In  this 
fomulation  the  requirement  of  corrmon  singular 
structure  corresponds  to  tho  condition  that  tho 
0 aid  0 have  the  same  normal  velocity  on  Sm, 
while  the  requirement  that  the  outer  boundary 
condition  for  0 shall  make  it  closely  equal  to 
0 on  S now  tokos  the  form 


Do  » B 


J“DJ  • 
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Once  Om  outer  boundary  condition  lia«  l>ocn 
established  wo  can  either  use  it  directly  to 
solve  the  coei|>l tented  tr.inslormed  version  of  the 
transonic  differential  pquation,  or  trnnsrurm  it 
back  to  tlw  orlrin.il  geometry  and  there  apply  it 
to  the  simpler  original  equation.  Tlie  choice  is 
ono  of  convenience  in  the  numerical  work.  In 
the  following  we  shall  l'  concerned  only  with 
establishing  the  outer  boundary  condition  in  the 
circular  goose  try. 

6.  Asymptotic  Approximation  for  Narrow  Slots 

The  fact  that  in  practice  slots  arc  usually 
narrow  as  compared  to  distances  between  slots 
will  be  the  basis  for  the  continued  analysis. 

In  this  situation  asymptotic  expansion  with  re- 
spect to  the  slot  width  as  a small  parameter 
suggests  itself  as  a useful  method  to  be  adopted. 
Obviously,  each  slot  will  need  its  own  'inner' 
expansion,  sealed  with  the  slot  width.  Tho 
•outor'  expansion  will  be  concerned  with  the 
overall  flow  on  the  scale  of  the  test  section 
radius. 


m 7 1-2r  coa(8-0*1*)  ♦ r*"  . 
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Fig.  3.  Inner  and  outer  flows. 

In  an  inner  expansion  the  slot  will  be  lo- 
cated alone  in  an  infinite  plane  wall  (Fig.  3). 

In  the  outer  expansion  the  test  section  wall 
will  be  solid  with  sinks  or  sources  located  at 
the  points  into  which  the  slots  have  contracted. 
The  total  flux  of  each  sink  or  source  will  of 
course  bo  equal  to  twice  the  corresponding  flux 
q (i  *=  1,2,... ,N)  through  the  slot  (positive  if 
into  the  plenum  chamber;  flux  unit  = U times 
test  section  radius).  These  fluxes  are  not  known 
In  advance  but  depend  on  tho  plenum  pressures  pID 
through  the  pressure  conditions  (5)  at  the  bound- 
aries S-g  and  Sp  . They  will  therefore  have 
to  be  obtained  by  matching  the  outer  expansion 
to  each  of  the  inner  expansions.  For  a survey 
of  matching  problems  involving  flow  through 
narrow  slots  seo  Ref.  12. 

Leaving  the  inner  expansions  for  later,  wo 
•hall  first  obtain  an  outer  representation  of  0, 
assuming  the  fluxes  q‘i’  to  bo  known.  This  will 
permit  us  to  analyse  0 and,  in  fact,  to  make 
a first  specification  of  tho  boundary  condition 
for  0. 

An  elementary  solution  with  a source  at  tho 
origin  and  a. sink  of  double  strength  on  the  unit 
circle  at  8 ■ 8'-^  (Fir,  <»)  is  given  by  ( yijn)  In  r- 
(l/njln#  , where  r'*'(r,8)  is  tlie  distance 
from  the  sink! 


'SOURCE 


Fig.  U.  A source  at  the  origin,  a sink  of 

double  strength  on  the  unit  circle. 

It  is  easily  shown  that  this  solution  lias  rero 
normal  derivative  on  the  unit  circle  (for  8$  8“’). 
We  can  therefore  immediately  write  down  tlie 
outer  representation  of  0 l 

0 = ± £ qli,lniii,e  G(r,0).  (12) 


The  undetermined  function  G is  harmonic  in  the 
unit  disk,  except  at  the  origin,  and  has  zero 
normal  derivative  at  tlie  outer  boundary.  Com- 
parison with  (9)  shows  that  A * £ q'**/2n  and 
that  ° ' 

G = D + £ (A.cosj9+  B .slnj9)(r-^+  r^).  (13) 

J=1  J J 

To  obtain  0 we  must  expand  the  rest  of  0 . 
It  is  easy  to  show  that 

In  i*1’.  - £ -j  cos[j(8-e<l))]rJ+  Otr^1),  (l«l) 


5 = ^(Sq<l))  • In  r + D ♦ 


♦ l£q(i,{  £-1  cos[j(8-ef°)]rJ}  + 

♦ £ (A.cosjS*  B,sinJ8)(r~^  + r^)  * 

J=1  J J 

- 0+4  " r<1)  * £ 4cos[j(e-e<1))]rJ}  . (15) 

% "i«i  *■  j=iJ  * 


This  gives  imncdiateiy  an  important  part  of 
tlie  attempted  outer  boundary  condition  for  0 1 

rmii  i *4  E *,U>  E cos[j(e-e(1,))V  . (16) 

"l«l  J-t  * 

This  expression  for  tlie  normal  derivative  of  0 
at  the  outer  boundary  docs  not  contain  explicit- 
ly the  coefficients  Aj  and  B.  , which  repre- 
sent tho  shape  of  the  model,  and  so  it  seems  to 
be  local  ns  well  as  homogeneous  in  tho  previous- 
ly defined  sense.  Ono  might  suspect  that  the  un- 
known fluxes  q*  , when  resolved,  will  bring  tho 
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modol  bark  into  tl*o  picture,  hut  since  they  ore 
to  be  obtained  by  matching  at  the  slots  we  ex- 
pect thorn,  or  rourso,  to  be  locally  determined 
in  terns  of  f , 

In  preparation  for  the  matching  tm  deduct 
from  (15)  the  following  iiuier  representation  of 
the  outer  representation  of  f at  the  slot 
point  k 1 

(k)  N 

lm1  (17) 


H W.-l  Inr'^-i^corf.e'11)),^,  • 

J=1J 

H(k,k)  = _l  £l  (ifltBl N)  , (18) 

J=1J 

with  r^,k*  denoting  the  distance  between  the 
slot  points  i and  k . The  dominating  part  of 
the  remainder  corresponds  to  a tangential  flow 
past  the  slot  point. 

# 7.  Inner  Representations  and  Matching 

The  inner  flows  might  be  analysed  either  in 
tho  transformed  or  the  original  cross-flow  plane 
since  tho  matching  will  involve  only  0 and  the 
q(i* , which  are  all  invariant  tinder  the  conformal 
snipping.  We  choose  to  work  in  the  original  plane 
as  being  clearly  the  simplest  alternative. 

As  explained  in  general  terms  in  Section  2, 
we  shall  work  with  a simplified  slot  flow  model 
(Fig.  l)  in  which  fast  air  from  the  test  section 
enters  the  plenum  chamber  as  a Jet  while,  in  re- 
verse flow,  the  fast  air  in  the  slot  returns  to 
tho  test  section  without  separation,  letting  in 
quiescent  air  behind  it  to  form  a longitudinal 
separation  bubble  at  plenum  pressure  Pp*4*  . 

Since  the  boundary  S between  the  two  masses 
of  air  is  not  precisely  defined  in  the  model,  we 
can  use  this  freedom  to  simplify  the  computation 
of  the  inner  representation.  More  specifically, 
we  siiall  assume,  at  least  in  the  steady  case, 
that  for  any  given  family  of  similar  slot  cross- 
sections  we  can  restrict  ourselves  to  one  basic 
cross-flow  solution,  taking  the  different  posi- 
tions of  one  of  its  material  curves  to  represent 
possible  free-boundary  curves  Sp  . Then  tho 
slot  flow  model  is  essentially  complete  as  soon 
as  we  prescribe  the  curve  at  which  the 

pressure  condition  is  to  be  satisfied  when  there 
is  a Jet  into  the  plenum  chamber  and  at  which  the 
Jot  is  to  be  tnken  to  sp]  i t when  tho  flow  revers- 
es. Introducing  Cartesian  coordinates  (z,y) 

(We-  5),  we  take  as  the  parameter  for  any  curve 
S_  its  intersection  y with  tho  y-axis;  y , 
the  maximum  valuo  of  yj,  , corresponds  to  SpoP. 

Now  lot  a (x)  bo  a measure  of  the  width  of 
the  slot.  The  velocity  level  in  the  slot  mid  its 
neighbourhood  is  given  by  »fk,/a  • With  this 
should  be  compared  the  velocity  set  up  by  a 
stream vise  variation  of  the  slot  width,  as  ex- 
pressed by  tho  normal  velocity  H in  (3).  As 
this  is  proportional  to  da/dx  , it  is  evidently 
two  orders  of  magnitude  smaller  in  the  limit 


Fig.  5.  Slot  geometry. 


a — 0 and  we  can  therefore  safely  neglect  it 
when  computing  the  cross  flow.  In  addition  to  a 
sink-like  flow  into  the  slot,  the  geometry  allows 
a flow  essentially  parallel  to  the  wall.  Its 
velocity  is  determined  by  the  outer  flow  and  is 
consequently  0(q(k*).  We  can  neglect  this  flow 
component,  too,  when  writing  down  the  inner  re- 
presentation of  0 to  the  lowest  order;  it 
should  bo  included  in  the  next  order  approxima- 
tion, however. 

We  first  solve  the  sink-flow  problem  with 
act,  q*W)  - | and  with  vanishing  potential  at 
y b z = 0.  (How  this  can  be  done  is  demonstrat- 
ed in  Appendix  2 for  the  particular  geometry  of 
5*)  hot  this  normalized  potential  be 
Q(z,y).  Then  our  inner  representation  is 

0 = *,k,+  q,k,-Q(z/a,  y/a)  , ( 19) 

wher  ' (x,t)  is  the  value  of  0 for  slot  k 

at  y = 2 = 0 (if  necessary  extrapolated  from 
Sp  by  Q ) . 

Far  away  into  the  test  section  Q has  the 
representation 


Q(z.y)  = - ^ in  rW  + R + °(^}k))  * 
re 

R = lim  («  t ; m rW  ) . 


The  corresponding  representation  of  0 is 

0 = 0^-  tn(r^/a)  * q^R  + 0(-^ ) , 

and  this  can  be  matched  to  the  outer  representa- 
tion (17)*  lienee 

% 

0,k'=  0 ( 1 ,6,k*  ;x,t)+  t qW H*i,k)  - 
iwl 


in  a + R)  . 


This  important  result  permits  us  to  relate  0 
locally  to  the  plenum  pressures. 

8.  The  Tlenum  Pressure  Condition  in  Steady  Flow 

So  far  the  analysis  applies  in  its  main  fea- 
tures to  unsteady  as  well  as  steady  flow.  We 
now  restrict  ourselves  to  tho  steady  case  (to 
return  later  for  a brief  discussion  of  the  un- 
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•toady  ease).  Ve  can  tlien  apply  tlie  pressuro 
condition  (5),  which  gives  at  r.  * 0,  y = yp  o 


where  v,,(y  /a)  Is  tho  value  of  tho  y-derivn- 
tive  of  the 'normal ized  potential  Q at  z = O on 
S„.  It  is  sufficient  to  satisfy  the  condition 
at  z « 0 if  tho  variation  of  flow  velocity  on 
the  boundary  curve  is  small,  ns  we  shall  assume 
in  concordance  with  our  simplified  description 
of  Sp  . In  order  to  keep  truck  of  y^fx)  when 
.it  is  smaller  than  yp0  one  has  to  solve  tho 
differential  equation 


- V Vp(y/R)  • (23> 

marching  downstream  from  x = xD , where  y = 0, 
and  taking  a new  initial  value  >’p0  each  time 
the  jot  splits  on  flow  reversal. 

Inserting  0 from  (19)  into  (22),  with  0^ 
from  (21),  we  obtain  tho  pressure  condition  in 

tho  form 

' & !«.. <H,x>  ."JW-M. 

1=1 

+ qW[ti(0  Yp/a)-  R - £ In  aj  ♦ 

♦ !(qWv/a)S  *-  6lk)  (k  = 1 N)  . (2!t) 


If  all  tho  fluxes  ql  (x)  are  known  this  is  an 
ordinary  differential  equation  for  $(1,8*^  j x ) 
which  can  be  integrated  together  with  (23), 
starting  with  a known  value  for  0 at  the  be- 
ginning of  the  slot. 

9.  On  Outer  Boundary  Conditions  for 
Transonic  Flows 

At  this  stage,  before  we  attempt  to  construct 
an  outer  boundary  condition  from  the  raw  material 
gathered  in  tho  preceding  sections,  it  might  be 
useful  to  consider  wrliat  sort  of  a boundary  con- 
dition we  should  like  to  have. 

Ve  are  now  concerned  with  steady  flows  only, 
•o  tho  differential  equation  with  wltich  to  use 
the  boundary  condition  is 

Af  - [M»  . (25) 

Typically,  with  those  unbounded  transonic  flows 
for  which  the  small-perturliation  analysis  is 
valid,  tho  right-liand  member  is  small  compared 
to  the  cross-flow  derivatives  constituting  the 
left-hand  member.  Ve  expect  this  to  be  true 
also  in  the  present  case,  at  least  to  the  extent 
that  «t  have  been  successful  in  designing  a slot- 
ted test  section  with  small  interference.  Tills 
means  of  course  tliat  the  cross  flow  in  any  plane 
X ■ constant  is  approx imately  volume  conserving 
so  that  as  much  volume  flux  as  entering  at  tho 
inner  boundary,  there  determined  by  Iho  shape  of 
the  modol,  must  leave  through  the  outer  boundary 
at  about  tho  same  cross  section.  Tills  ocnsutulos  n 


strong  restriction  on  <fy,  nt  the  outer  boundary 
in  terms  of  <f(|  at  the  inner  boundary,  suggest- 
ing that  we  might  run  into  difficulties  if  we 
try  to  prescribe  tho  normal  derivative  at  the 
outer  boundary. 

These  considerations  come  sharply  Into  focus 
in  tlic  case  of  axi  symmetric  flow  in  tho  slonder- 
body  approximation,  that  is  with  the  riglit-hand 
member  of  (25)  neglected.  Then  the  perturbation 
potential  <p(x,r)  takes  the  form 

rVr  - »'(*).  f = s'(x)  • In  r ♦ f(x,l)  , 

whore  s(x)  is  essentially  the  cross-sectional 
.area  of  the  body.  In  this  case  the  radial  vel- 
ocity is  everywhere  exactly  determined  by  the 
inner  boundary  condition,  at  the  body,  and  only 
the  potential  itself  can  bo  prescribed  at  an 
outer  boundary  (at  r=  1,  say).  Therefore  we  ob- 
tain a reasonable  slender-liody  solution  if  the 
outer  boiuidaty  is  a jet  boundary  at  prescribed 
pressure,  or  a ventilated  wall  which  responds  to 
the  radial  velocity  s'  (x)  with  a well-defined 
pressure  distribution.  On  the  other  liand,  with 
a solid  wall , at  which  zero  normal  derivative 
must  be  prescribed,  there  is  110  solution  at  all. 
Then  the  riglil'-hand  member  must  be  retained  in 
(25),  accounting  for  tho  raised  blockage  inter- 
ference level  typical  of  choking. 

Ve  shall  therefore  attempt  to  obtain  our 
outer  boundary  condition  in  the  form 


on  S„  1 


? = ? tfn]  . 


whore  jf"  [ ) is  a regular  functi  onal  over  Sw 
(it  becomes  singular  in  the  limit  of  a solid 
wall,  of  course).  This  choice  carries  with  it 
tho  suggestion  that  in  solving  (25)  by  an  itera- 
tive numerical  method  one  should  march  back  and 
forth  between  the  model  and  the  outer  boundary, 
arriving  each  time  at  the  outer  boundary  with  an 
improved  normal  velocity  to  give,  by  (26),  an  im- 
proved pressure  distribution  to  be  carried  back 
to  the  model.  It  is  reassuring  to  be  able  to 
note  that  employing  such  a scheme  speeds  up  con- 
vergence, ill  particular  in  subsonic  regions.13 
Our  philosophy  also  suggests  that  the  finite- 
difference  approximation  f or  the  left-hand  member 
of  (25)  should  be  a conservative  one. 

Tho  choice  of  (26)  should  not,  however,  pre- 
clude interest  in  the  inverse  functional,  jt'1  , 
giving  !pn  ill  terms  of  ip.  Actually,  if  ve  are 
to,  adopt  the  idea  of  a 'self-correcting'  test 
section,  presently  being  pux-sued  in  several  lab- 
oratories, then  a viable  scheme  might  bo  to  meas- 
ure the  wall  pressure  distribution  along  slats, 
essentially  <p , and  hence  compute  the  normal  vel- 
ocity tpn  rather  than  trying  to  measure  it.  For 
this  $~x  would  be  needed  (together  with  a pro- 
cedure for  obtaining  (p  from  <p  f see  Sec- 
tion 12).  One  would  next  liave  to  compute  the 
unbounded  flow  outside  Sw,  using  ifL  In  an 
timer  boundary  condition  and  obtaining  an  esti- 
mate of  what  ip  should  be  on  Sw  in  order  to 
be  free  of  interference.  Using  again  the  same 
qi,,  , onn  would  then  adjust  the  test  parameters 
so  as  to  make  ip,  computed  hy  (26)  with  tho  ad- 
justed 7 , equal  to  ip  obtained  from  tho  oulor 
flow.  Repeating  this  cycle  ns  long  os  tho  re- 
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qulrcd  wall  nrtjustmi'iils  can  bo  realized,  employ- 
ing alternately  > 1 ui«l  J , one  .nlioulil  end  up 
with  * slot  toil  wall  oT  minimum  interference. 
Running  the  cycle  Uiekwurtls,  adjusting  ^i)l  for 
equality  rather  tlion  if  , would  seem  to  l»e  In- 
viting difficulties. 

10.  Constructing  the  Outer  Boundary  Condition 

According  to  the  prescription  of  Section  k 
for  tlie  boundary  condition  on  Sw  , ij>  is  to  be 
related  to  qy,  in  the  same  way  as  j is  re- 
lated to  $n . To  the  extent  that  we  have  estab- 
lished this  relationship  it  is  contained  in  Equa- 
tions  ( 16) , (23)  and  ( 2*» ) . 

We  take  ( 16)  as  a starting  point  since  it 
contains  the  normal  derivative,  the  ar/jinnonl. 
function  of  O'  . Substituting  <j>(x,r,9)  for 
$(r,9jx)  and  specializing  9 to  the  slot  posi- 
tions, vc  obtain 

t cos[j(9lkLe,t))]}c,(iJ 

i=1v  J=1  J 

“ "•?,.(*.  (K  = 1,2 N)  . (27) 

Wiis  a linear  system  which  determines,  for  each 
x,  the  slot  fluxes  in  terms  of  tfr  • The  coef- 
ficient matrix  is  independent  of  x and  can  be 
inverted  numerically  and  stored  as  soon  as  the 
■lot  locations  have  been  decided.  Tims  it  is 
very  simple  to  calculate  the  slot  fluxes  from  <fr . 

Next  vc  integrate  the  equations  (23)  to  de- 
termine the  flow  penetration  depth  in  the  slots. 
Tills  is  straightforward,  taking  one  slot  at  a 
time.  Now  tho  slot  geometry  is  involved,  how- 
ever, so  the  integration  must  be  repeated  each 
time  we  make  a wall  adjustment. 


*x(x„,..9,k,)..[i(vpllm  a 
o 


£*)\  6W] 


to  bo  wel  1-dolermined  since  llm(q  /a)  will 
nover  be  known  wi th  any  precision.  This  gives 
us,  perhaps,  the  opportunity  to  specify 
limfqM/a) , which  is  also  needed  in  Eq.  (23), 
in  such  a way  that  takes  the  value  corre- 

sponding to  the  chosen  entrance  Mach  number  M , 
thus  avoiding  a spurious  velocity  jump  at  x = x . 
in  practice  the  achievement  of  a smooth  ontranc! 
flow  might  require  a proper  choice  of  plenum 
pressures. 


In  order  to  complete  tho  construction  of  tho 
functional  3"  wc  have  only  to  define  if  between 
slots.  In  accordance  with  our  method  of  analysis 
this  is  done,  in  each  plane  x = constant,  by  tri-  1 
gonometric  interpolation  of  order  u . Obviously  I 
u must  be  chosen  so  as  to  make  2o  + 1 equal  to 
N at  most  ; usually  it  will  be  smaller.  Tills 
means  that  usual ly  wc  cannot  satisfy  the  condi- 
tions (28)  precisely.  A least-squares  fit  could  v 
then  be  used,  which  again  calls  for  storing  a 
pre-conijmted  constant  matrix.  The  deviations,  of 
course,  correspond  to  higher  order  components  to 
be  filtered  out.  Since  both  (.23)  and  (28)  are 
non-linear  with  respect  to  q(J)  such  components 
will  arise  even  if  not  present  in  the  tfC.  used 
in  (27).  ^ 


llic  inverse  functional  is  not  so  easily  con- 
structed, due  to  the  non— lineari ty  of  the  equa- 
tions from  which  are  determined  the  fluxes  q(*-J 
correspond iug  to  a given  potential  on  Sw  . An 
iterative  scheme  must  be  set  up.  Once  the 
fluxes  ore  known,  the  normal  derivative  is  im- 
mediately obtained  from  (27)  at  the  slot  posi- 
tions. The  construction  is  again  completed  by 
trigonometric  interpolation. 


Prom  (2*1),  finally,  we  obtain  the  value  of 
the  potential  at  the  slot  positions, 

v(x,1,eW)  = ,po_  [Q(0,yjj/a)-R  -i  In  a)q(k^- 

- _j  (j(  ,M  « 


(ic  — 1|.*.|N), 


(28) 


satisfying  at  x = xc  the  upstream  condition 
for  $ assumed  in  Section  3.  The  righl-liand 
member  is  now  completely  known.  It  would  re- 
main valid  if  we  were  to  allow  the  plenum  pres- 
sures to  vary  with  x . 

It  is  primarily  Equation  (28)  we  must  ana- 
lyse when  wo  want  to  adjust  slot  flow  parameters 
so  as  to  produce  the  wall  pressure  distribution 
corresponding  to  interrei-eiico-froo  flow.  There 
will  always  bo  interference  around  the  entrance 
section  at  x = x0  . The  test  parameters,  inclu- 
ding tho  direction  of  the  unbounded  flow  wo  are 
trying  to  simulate,  should  be  selected  in  such 
• way  as  to  make  this  interference  small  and 
Ideal.  Also,  note  the  singular  behaviour  of  $ 
•It.  *„»  vljore  not  only  n vanishes  but  also 
q*KJ  (sinco  Grad  9 = 0 over  the  whole  plane 
*■*.)•  Ono  cannot  expect 


11.  Symmc  trie  Flow 

As  a simple  application  consider  the  case  of 
an  axisyrxnetric  body  along  tho  axis  of  a circular 
test  section  with  N uniformly  distributed  iden- 
tical slots  (N  > 1).  From  symmetry,  the  flow  is 
periodic  with  respect  to  9 with  period  2u/N  . 
The  potential  q>  will  therefore  be  independent 
of  8 to  order  N-1  in  a trigonometric  expan- 
sion. Since  N-1  is  larger  than  (n-1)/2,  tho 
upper  hound  for  u , we  conclude  that  <p  is 
independent  of  8 . Actually,  since  the  flux 
through  each  slot  is  given  by 


Q=*d  ip_(x,l),  d = “ (=  arc  length  between  slots) 

N (29) 

and  is  independent  of  u , we  will  get  tho  same 
ip(x,r)  for  all  permitted  values  of  U< 


Equation  (29)  takes  the  place  of  (27).  It 
remains  to  consider  (23)  and  (28).  Noting  that 
in  the  present  case 


N 

£ II 
i=1 


(i,k) 


- InN 


■wo  obtain 


(30) 


r. 


/ 
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and  x 

9(x,l)  > f0- Kvr(*.')-  ][^(v  ^ Tr)  (3') 

X p 

with  ° 

K “ d[n  ™ " n t,^_R+  Q(°*y'n)]  * (33) 

E<i>iaiion  (31)  nprcos  with  the  classical  for- 
mula, augmented  by  u quadratic  cross-flow  term 
in  the  manner  of  W.W.  Wood14  • It  is  shown  in 
Appendix  2 that  Equation  (32)  is  a second  order 
approximation  for  small  a/d  to  the  classical 
ns  well  as  the  general Azed 4 K. 

It  may  seem  remarkable  that  9 is  independ- 
ent of  0 , It  does  not  ii»'*an,  however,  that  the 
vail  interference  is  totally  independent  of  the 
number  of  slots,  not  oven  for  a fixed  9 . In 
adding  slots  wo  suppross  successively  hi/^her 
harmonics  which  were  neglected  in  $ (but  pres- 
ent in  9 ) when  the  slots  were  fewer, 

Tlie  outer  boundary  condition  9 = as 

constituted  by  Equations  (30),  (30  and  (3-)  has 

been  tested  quite  extensively  4 with  the  numeri- 
cal method  of  Kef,  13,  No  difficulties  were 
found  except  vhc»re  the  slot  width  was  very  small, 
smaller  than  normally  used,  in  which  case  the 
Convergence  l»cca/ne  slow.  This  is  of  course  not 
surprising  since  S?  is  singular  in  the  limit  of 
a solid  wall . 

12.  Analysis  01"  Vail  Pressures 

For  a further  appl  ical. i on,  assume  that  wo 
want  to  determine  9 on  Sv  by  measuring  the 
pressure  distribution  along  slats,  *at  constant 
0 « 0W  t 0d)  say.  We  first  integrate  (*l)  along 
the  slat  to  get  9 (the  quadratic  cross-flow 
term  can  he  neglected  in  this  case).  Ihen  from 
Equations  (7)  and  (15)  wo  have 

9 = 9 ^ 0 — 0 

N u 

= rW.*  -1  cos [j (ev-e<y)l } (33) 

wliore  rw^  is  the  cross-flow  distance  from  a 
pressure  tap  to  the  slot  point  i.  In  the  case 
of  symmetric  flow  with  9W  half  way  between 
slots  tliis  reduces  to  the  result  obtained  in 
Ref.  *s. 

The  application  of  (33)  is  straightforward 
if  we  know  the  fluxes,  but  we  don't  in  the  inter- 
esting case  when  we  want  to  detennino  q(*)  arxl 

9p  by  applying  to  9 , However,  0-0 

has  the  cliaracter  of  a small  correction,  so 
there  is  little  doubt  that  an  iterative  scheme, 
using  successively  improved  estimates  for  , 

will  converge  rapidly.  This  is  not  really  very 
much  of  a complication.  $’x  has  to  be  computed 
iteratively  anyway. 

• 13.  Unsteady  Flows 

In  the  first  seven  sections  the  analysis 
applies  in^its  main  features  also  lo  unsteady 
flows.  In  the  pressure  condition  (5)*  however, 
wo  must  add  9^/0  to  9X  , recognizing  at  the 


some  time  that  the  free  surface  Sp  might  move 
and  that  the  plenum  pressures  6'  ' are  no  longer 
known  in  advance  since  pressure!  waves  may  propa- 
gate into  H10  plenum  cliamhor,  Tlio  simplified 
picture  of  the  slot  flow  introduced  in  Section  7 
should  be  extended  to  include  a potential  which 
describes  the  cross  flow  due  to  the  motion  of 
Sp  , in  addition  to  the  potential  Q which  de- 
scribes the  cross  flow  due  to  the  obliqueness  of 
the  flow.  (This  approach  means,  of  course,  that 
we  are  assuming  the  slot  width  to  be  small  com- 
pared to  typical  wave  lengths  of  the  unsteady 
flow.)  The  additional  cross  flow  is  somewhat 
different  in  tliat  it  involves  the  quiescent  air 
of  the  plenum  chamber  as  well  as  the  fast  air  in 
the  slot.  Tliis  will  also  add  a new  contribution 
to  the  right-hand  member  of  Eq.  (23)  for  Sp  , 
the  left-hand  member  of  which,  of  course,  now 
becomes  [d/dx  *■  (l/u)d/dt]>'p  . All  these  effects 
combine  to  make  the  pressure  condition  corros|ionci- 
ing  to  (2^l)  rather  complicated  in  its  final  form. 

In  order  to  account  for  fluctuations  in 
we  must  analyse  the  wave  propagation  in  the  plen- 
um chamber,  a forbidding  task.  The  simplest  case 
arises  if  the  plenum  chamber  is  so  large  that  only 
waves  propagating  from  the  slots  into  the  plenum 
chamber  contribute  to  the  pressure  fluctuation 
on  Sj)Q  + Sp  • Those  waives  correspond  to  the  ad- 
ditional, non-steady  part  of  the  slot  flow  and 
so  can  be  analysed,  assuming  a known  steady  pres- 
sure deeper  inside  tlio  plenum  chamber.  In  the 
next  more  comp]  icat ed  c.ase  one  would  have  to 
account  for  varying  plenum  pressure  in  the  inaiiier 
of  the  classical  Helmholtz  resonance  analysis. 

It  scorns  that  in  these  two  cases  there  is  a rea- 
sonable chance  to  be  able  to  complete  tlio  ana- 
lysis and  construct  a homogeneous  boundary  con- 
dition for  the  test  section  flow.  It  will  per- 
haps bo  necessary  to  linearize  the  problem, 
assuming  the  unsteadiness  to  be  a small  perturba- 
tion of  a steady  flow. 

1*1 , Wall  In terference  Corrections. 

Test  Sec tion  Calibration. 

Suppose  we  have  succeeded  in  adjusting  the 
test  parameters  so  as  to  make  the  interference 
negligible  at  the  model.  Then  tlirce  wall  correc- 
tions of  classical  typo  are  ininediatel  y available 
from  this  process:  a Mach  number  correction 
(=  M-M0),  an  anglc-of-.attack  correction  (=  the 
ancle  of  attack  of  the  test  section  with  respect 
to  the  direction  of  the  unbounded  flow  simulated), 
and  an  anglc-of-yaw  correction  (=  the  correspond- 
ing angle  of  yaw  of  the  test  section). 

« The  other  classical  corrections,  those  for 
induced  buoyancy  and  flow  curvature,  are  rightly 
absent.  As  soon  as  they  arc  needed  there  is  a 
distortion  of  the  pressure  distribution  over  the 
model  tliat  cannot  bo  tolerated.  In  contrast, 
the  former  corrections  are  not  associated  with 
any  such  distortion  and  should  be  jiennittod  to 
be  large  if  it  helps  in  reducing  the  interfer- 
ence. Tliis  is  "the  principle  of  minimizing 
interference  rather  than  corrections"15  , iripl  ying 
the  concept  of  a "correctable-interference  trans- 
sonic wind  tunnel"18. 

There  might  bo  other  errors  present  in  the 
tost  section  flow,  errors  which  are  not  accounted 
for  when  computing  the  wall  interference.  Dis- 
turbances from  tlio  entrance  section  upstream  of 
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tlm  slots,  oxiai  gradient «»  set  up  by  tho  wall 
boundary  layers  or  by  Improper  alignment  of  tho 
walls,  or  upstream  influence  from  tho  sting 
arrangement  aro  examples.  The  main  tool  for 
handling  such  errors  is  "calibration”,  running 
tho  tost  section  empty  with  the  hulls  atxl  plenum 
pressures  set  for  uniform  flow.  Such  errors 
should  be  eliminated  by  proper  adjustment  of  the 
teat  section,  not  by  introducing  corrections  to 
tho  tost  data.  The  practice  of  using  the  Mach 
number  of  the  empty  test  section,  calibrated 
against  the  plenum  pressure,  as  a free- stream 
reference  Mach  number  for  model  tests  does  not 
seem  to  make  sense  in  the  context  of  the  present 
theory, 

15.  Concluding  Remarks 

The  present  inviscid  theory  of  wall  inter- 
ference in  slotted  test  sections  generalizes  to 
tliree-dimensions  the  theory  for  two-dimensional 
tests  developed  on  a classical  basis5’  14  in 
Ref.  k.  In  analysing  the  local  flow  at  each 
slot  separately,  tho  structure  of  the  theory  is 
such  as  to  facilitate  later  inclusion  of  correc- 
tions to  account  for  viscous  effects  inside  the 
slots  and  the  plenum  chamber.  A first  attempt 
in  this  direction  has  already  been  made4  and 
V more  extensive  comparisons  with  experiment 
are  completed  it  might  become  possible  to  extend 
and  delineate  the  area  where  the  inviscid  theory, 
with  or  without  corrections,  might  be  used  with 
confidence.  The  interaction  of  shock  waves  with 
the  slot  flow  must  bo  studied,  in  particular. 

Meanwhile  the  inviscid  theory  enn  be  used 
for  running  numerical  experiments.  These  will 
show  how  accurately  one  must,  describe  the  action 
of  the  slotted  wall  in  different  types  of  appli- 
cation, and  what  wall  adjustment  facilities  one 
must  provide  in  order  to  eliminate  tho  wall  in- 
terference. They  will  also  help  developing 
strategies  for  efficient  use  of  adjustable  slot- 
ted walls  in  future  wind  tunnels. 
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Appendix  1 — Symbols 

a(x)  Slot  width}  Fig.  5 

d Arc  length  between  slots;  Eq.  (29) 

Grqd.^  Differential  operators  in  plane  x = 
constant 

3",  3-  Functionals  on  Sw  ; Eq.  (26) 

Eq.  (18) 

K(x)  Coefficient  in  slotted-wall  boundary 

condition;  Eqs.  (32),  (A-8)  and  (A-12) 

M Mach  number  of  reference  flow 

Me  Entrance  Mach  number;  Fig.  2 

N Number  of  slots 


p Pressure 


, Prossure  of  flow  at  Mach  number  M 

y 

r* 
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- 


f' 


Q(*»y)  Normalized  slot  flow  potential;  Eq.  ( 19) 
qllU*,t)  Flux  through  slot  1 ; Eq.  (12) 

R Eq.  (20) 


rW 

jJi.k) 


t 

U 

v (y  ) 

p'  v' 


m 

y.« 

yp 


1 

6W 


<p»? 


Radius  vector  in  cross-flow  plane 

Polar  coordinates  in  cross-flow  plane 

Distance  from  slot  point  i ; Fig.  h, 

Eq.  (11) 

Distance  lietwecn  slot  points  i and  kt 
Eq.  (18) 

Cylindrical  surface  enclosing  the  model 
and  its  wake;  Fig.  2 

Free  siu'facc  be  tween  fast  air  and 
quiescent  plenum  air;  Fig.  1 

Surface  at  plenum  pressure  across  jet 
on  plenum  side  of  slot;  Fig.  1 

Outer  boundary  of  test  section;  Fig.  2 


Time 

Flow  velocity  at  Mach  number  M * 

= &q/*y  at  y = yp,  r.  x=  0;  Eq.  (23) 
Distance  along  tunnel  axis;  Fig.  2 
location  of  entrance  section;  Fig.  2 
Cartesian  coordinates  at  slot;  Fig.  5 

Coordinate  of 
Fig.  5 

Coordinate  of 
Fig.  5 

Ratio  of  specific  heats 

Non  dimensional  plenum  pressure  dif- 
ference; Eq.  (5) 

Perturbation  velocity  potentials; 

Eqs.  (l)  and  (6) 

Harmonic  functions  in  cross-flow  plane; 
Eq.  (7) 


Sp  on  slot  centre  line; 


S on  slot  centre  line; 


W Highest  order  considered  in  harmonic 

analysis  of  interference;  Eq.  (8) 

Density  of  flow  at  Mach  number  M 


Appendix  2 

Simplified  Analysis  of  Slot  Flow 


As  a basis  for  analysing  the  slot  flow  for 
the  geometry  of  Fig.  5«  consider  the  symmetrical 
flow  or  unit  flux  from  a half  plane  into  a slot 
of  unit  width  and  unlimited  depth  (sec  sketch). 


Ve  shall  determine  the  velocity  potential  Q(y,z) 
normalized  to  vanish  at  tlio  origin. 

Tills  is  accumpl  1 shed  l>y  mapping  the  flow 
region  conformally  onto  n half  plane  (-tr  * s o 
in  sketch),  transforming  the  slot  extremity 


(Z  = ioc)  into  tlio  origin  and  leaving  the  far 
field  undisturbed.  Die  required  complex  trans- 
formation is 


♦ iy 


."/s’— 0»  i O In  ^:2t7tlco^ 

* s • 


0 - V • (A-0 

wlierc  (»,$)  are  polar  coordinates  in  tile  trans- 
formed piano  C.  Tlio  exterior  wali  and  the  slot 
walls  are  mapped  onto  the  real,  axis  with  the 
comer  points  at  ±0  = ± l/n  . The'  origin  goes 
into  a point  on  the  negative  imaginary  axis,  at 
distance  s from  the  origin,  say.  The  fiow  is 
that  of  a sink  at  the  origin,  hence 


Q = 


2 

n 


(A-2) 


71iis  gives,  according  to  Eq.  (20),  R = (l/n)lnso< 

As  the  material  ciu-ve  to  represent,  in  its 
different  positions,  possible  free— boundary  curves 
sm  • vo  choose  one  which  far  into  the  slot  is  a 
straight  line  across  the  slot.  In  the  trans- 
formed plane  the  curve,  when  close  to  the  origin, 
is  a circle.  In  order  to  determine  its  develop- 
ment ve  can  integrate  along  streamlines  (rays)  in 
the  transformed  plane.  Points  on  the  same  curve 
must  have  the  same  value  of  the  parameter 

s 

*(».*)  = J ld|l  s ds  » (A-3) 

c 


where  c is  a small  positive  number.  After  the 
curves  have  been  determined  in  the  C-plnne  they 
have  lo  be  transformed  back  into  the  Z-planc. 
Ihis  has  not  been  carrier!  out  yet. 


For  our  immediate  purpose  it  is  sufficient 
to  consider  the  intersection  y of  Sm  with 
the  axis  of  symmetry.  This  point  is  mapped  on 
the  negative  imaginary  axis,  at  the  distance  s 
from  the  origin,  say.  The  value  of  sD  is  ob-* 
tained  from  the  equation  P 


+ o*t  a In 


(A-'«) 


*hich  gives,  in  particular,  s_  ■ s0  • 0.2110  for 
Yp  = 0.  Once  -Mp(yp)  is  known,  it  Is  strajglit- 
forvanl  to  compute  Q and  vp  1 


11 


„ . , i . d In  a 

9(0,y  )-R*»-  in  * , v * r.  (A-5) 

P n p’p  n«ly  * ' 

P 

If  yp  is  well  ins  i«'e  thn  slot  w»-  cmi  Coke 
?n/®  tx*  «tiall,  Ex)Kii^im/'  the  square  loot  in 
(A-q),  the  following  accoihI  order  approximation 
la  obtained: 


2 -2{ny  ♦!) 

■ 1 - In  — «ay  40  p 
P n p 


(A-6) 


Hence 


Q(Otyp)  - R - 1 (I  . In  2)4  yp*l  .“2("V,)t 

(A-7) 


. , -2(ny  4 1 ) 

VV  - > - sc  p 


Inserting  into  Eq.  (32),  we  obtain 
K“d[n  in™  + £ 0“  |)+  Vp/a  +1  c 


(A-8) 

in  complete  agreement  (to  second  order  in  a/d  ) 
with  the  result  of  Ref  A,  at  least  for  yp/a>0.1. 

This  result  is  valid  if  S is  inside  the 
slot.  If  the  flow  is  leaving  the  slot  as  a jet 
into  the  plenum  chamber  tile  analysis  in  Ref.  h 
shows  that  Eq.  (a-8)  is  still  valid  if  y/a  is 
cl  tanged  into  i/a  4 0.22,  where  I is  the  depth 
of  the  slot.  In  the  present  terminology  this 
means  simply  that  we  shall  take 


When  there  is  a 'bubble' , for  y < O , Eq. 
(A-^l ) can  again  be  solved  approximately.  The 
re sul i , for  Sp  » o , is 


s m - y 
P TP 


' 8n*  y » ' ’ 


hence 


p 

P " V l»ti3  y3  ' 


(a- 10) 


(A- 11) 


This  is  consistent  with  S becoming  approxi- 
"lately  a semicircle  with  fcentre  at  the  origin 
as  soon  ns  ly.,|  >a.  The  corresponding  result 
for  K is 


K 


<[i 


In 


d 

2"lypl 


1 1 

8n*  (y/a)*  J 


(A— 12) 


Finally,  it  might  be  noted  Hint  the  preced- 
ing analysis,  with  very  minor  modifications, 
applies  also  to  a slot  located  in  a right-angled 
corner.  One  has  only  to  make  the  Mis  of  sym- 
metry into  a wax  a and  everything  else  follows. 


*po  | 

* - a + °- 


22 


(A-9) 
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